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ABSTRACT 

We present new computations of the equilibrium and non-equilibrium cooling efficiencies and ioniza- 
tion states for low-density radiatively cooling gas containing cosmic abundances of the elements H, He, 
C, N, O, Ne, Mg, Si, S, and Fe. We present results for gas temperatures between 10 4 and 10 8 K, assum- 
ing dust-free and optically thin conditions, and no external radiation. For non-equilibrium cooling we 
solve the coupled time-dependent ionization and energy loss equations for a radiating gas cooling from 
an initially hot, > 5 x 10 6 K, equilibrium state, down to 10 4 K. We present results for heavy element 
compositions ranging from 10~ 3 to 2 times the elemental abundances in the Sun. We consider gas cool- 
ing at constant density (isochoric) and at constant pressure (isobaric). We calculate the critical column 
densities and temperatures at which radiatively cooling clouds make the dynamical transition from 
isobaric to isochoric evolution. We construct ion ratio diagnostics for the temperature and metallicity 
in radiatively cooling gas. We provide numerical estimates for the maximal cloud column densities for 
which the gas remains optically thin to the cooling radiation. We present our computational results 
in convenient on-line figures and tables (see http://wise-obs.tau.ac.il/~orlyg/cooling/). 

Subject headings: ISM:general - atomic processes - plasmas - absorption lines - intergalactic medium 



1. INTRODUCTION 

The collisionally controlled ionization states and ra- 
diative cooling rates of hot (10 4 — 10 s K) low-density 
"coronal" gas clouds are crucial quantities in the study of 
the diffuse interstellar and intergalactic medium, includ- 
ing hot gas in supernova remnants, the Galactic halo, 
active galactic nuclei, galaxy groups and clusters, and 
warm/hot "baryonic reservoirs" in intergalactic clouds 
and filaments. In this paper we present new computa- 
tions of the equilibrium and non-equilibrium cooling effi- 
ciencies and ionization states for radiatively cooling gas 
clouds with temperatures between 10 4 and 10 8 K. Our 
main emphasis is on the non-equilibrium cooling and ion- 
ization states that occur when an initially hot gas cools 
radiatively below ~ 5 x 10 6 K. Below this temperature, 
cooling can become rapid compared to ion-electron re- 
combination, and the gas at any temperature tends to re- 
main "over-ionized" compared to gas in collisional ioniza- 
tion equilibrium (CIE). Because the energy losses below 
~ 10 7 K are dominated by atomic and ionic line emissions 
from many species, the non-equilibrium cooling rates also 
differ (are suppressed for over-ionized gas) compared to 
equilibrium cooling. The "recombination lag" and depar- 
tures from equilibrium depend on the abundances of the 
heavy elements, and we present results for compositions 
ranging from 1CP 3 to 2 times the elemental abundances 
in the Sun. We include H, He, C, N, O, Ne, Mg, Si, S, 
and Fe, and we calculate the ion fractions as functions of 
temperature, for both equilibrium and non-equilibrium 
cooling, for all ionization stages of these elements. 

The non-equilibrium radiative cooling of hot (> 10 4 
K), highly- ionized gas, with emphasis on the associ- 
ated time-dependent ionization of the metal ions, is a 
classical problem first investigated by Kafatos (1973), 
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and subsequently by several authors including Shapiro 
& Moore (1976), Edgar & Chevalier (1986), Schmut- 
zler & Tscharnuter (1993), Sutherland & Dopita (1993), 
and Smith et al. (1996). Time dependent cooling in 
colder (T < 10 4 K) neutral hydrogen gas was stud- 
ied earlier by Bottcher et al. (1970), Jura & Dalgarno 
(1972), and Schwarz et al. (1972) (see also Dalgarno & 
McCray [1972]). These investigations demonstrated that 
below a few 10 6 K recombination lags can indeed lead 
to significant departures from CIE states and the asso- 
ciated equilibrium cooling rates. Computations of hot 
gas cooling efficiencies and ionization states assuming 
CIE have been presented in many papers dating back to 
House (1964), Tucker & Gould (1966), Allen & Dupree 
(1969), Cox & Tucker (1969), Jordan (1969), followed by 
Raymond et al. (1976), Shull & van Steenberg (1982), 
Gaetz & Salpeter (1983), Arnaud & Rothenflug (1985), 
Boehringer & Hensler (1989), and more recently Suther- 
land & Dopita (1993), Landi & Landini (1999), and Ben- 
jamin et al. (2001). 

In this paper we reexamine the fundamental problem 
of non-equilibrium ionization in a time-dependent radia- 
tively cooling gas. We focus on the behavior for pure ra- 
diative cooling with no external sources of heat or pho- 
toionization. Recombination lags and non-equilibrium 
ionization also occur in an initially hot gas undergoing 
rapid expansion and adiabatic cooling (Breitschwerdt & 
Schmutzler 1994; 1999). We do not consider such dy- 
namical effects here. 

Our work is motivated by ultraviolet (UV) and X- 
ray metal absorption-line spectroscopy of hot gas in the 
Galactic halo (e.g., Sembach & Savage 1992; Spitzer 
1996), including, more recently, in discrete ionized "high- 
velocity clouds" (Savage et al. 2002; Foxet al. 2004, 2005, 
Fox, Savage, & Wakker 2006, see also Sternberg, Mc- 
Kee, & Wolfire 2002; Gnat & Sternberg 2004; Mailer & 
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Bullock 2004). We are also motivated by recent Hubble 
Space Telescope (HST), Far Ultraviolet Spectroscopic Ex- 
plorer (FUSE), Chandra, and XMM-Newton, detections 
and observations of hot gas around the Galaxy (Collins 
et al. 2005; Fang et al. 2006), and/or of gas that may be 
part of a 10 5 -10 7 K "warm/hot intergalactic medium" 
(WHIM) (Tripp, Savage, & Jenkins 2000; Shull, Tum- 
linson, & Giroux 2003; Richter et al. 2004; Sembach et 
al. 2004; Soltan, Freyberg, & Hasinger 2005; Nicastro et 
al. 2005; Savage et al. 2005, Williams et al. 2006; see how- 
ever Rasmussen et al. 2006). This includes the tempera- 
ture range where departures from equilibrium cooling arc 
expected to be important. An intergalactic shock heated 
WHIM is expected to be a major reservoir of baryons in 
the low redshift universe (Cen & Ostriker 1999; Dave et 
al. 2001; Furlanetto et al. 2005). Important ions for the 
detection of such gas include C IV, N V, O VI, O VII, 
O VIII, Si III, Si IV, Ne VIII, and Ne IX. 

The early discussions of Kafatos (1973) and Shapiro & 
Moore (1976) focused on time-dependent isochoric (con- 
stant density) cooling at a single metallicity. Suther- 
land & Dopita (1993) calculated isobaric (constant pres- 
sure) cooling efficiencies for a range of metallicities, 
but presented limited results for the time-dependent 
ion-fractions, and did not consider isochoric cooling. 
Schmutzler & Tscharnuter (1993) considered isochoric 
cooling for a range of metallicities, and presented ion- 
fractions only for solar metallicity gas. They presented 
limited results for isobaric cooling. In our calculations, 
we use up-to-date atomic data, and we present complete 
results for the time-dependent ion fractions and cooling 
efficiencies, for both isochoric and isobaric evolutions, 
for a wide range of metallicities. Because of the ii PdV'' 
work, isobaric cooling is less rapid than isochoric cooling, 
so departures from equilibrium are somewhat smaller in 
the isobaric case. 

We consider metallicities ranging from 10 -3 to twice 
the solar abundances of the heavy elements. In our 
calculations we assume that any initial dust is rapidly 
destroyed (e.g. by thermal sputtering) on a time-scale 
that is short compared to the initial cooling rate (Smith 
et al. 1996), and we consider dust-free cooling at fixed 
gas-phase elemental abundances. The sensitivity of the 
cooling efficiencies to the metallicity was studied by 
Boehringcr & Hcnsler (1989) for equilibrium cooling, 
and by Schmutzler & Tscharnuter (1993) and Suther- 
land and Dopita (1993) for non-equilibrium conditions 
as well. Time-dependent effects are enhanced (dimin- 
ished) as the cooling efficiencies are increased (reduced) 
for larger (smaller) gas-phase abundances of the heavy 
elements. Thus, recombination lags are more significant 
as the metallicity is increased. For a low-density gas, the 
non-equilibrium ion fractions and cooling efficiencies are 
independent of the gas density or pressure, because the 
ratio of the cooling and recombination times is indepen- 
dent of density. The density scales out of the problem. 
If the gas is assumed to cool, at either constant density 
or constant pressure, from an initial equilibrium state at 
high temperature, the only remaining free parameter in 
the problem is the metallicity Z . 

We present results for the optically thin limit in which 
rcabsorption of diffuse "self-radiation" may be neglected. 
We verify previous findings that trapping of hydrogen 
and helium recombination radiation, in optically thick 



clouds has only a small effect on the cooling and result- 
ing ionization states. We also provide estimates of the 
maximal cloud column densities for which reabsorption 
of radiation from the metals (lines and continua) may 
be neglected in computations of the gas cooling efficien- 
cies. Elsewhere, we will consider the effects of exter- 
nal fields on the non-equilibrium evolution of radiatively 
cooling gas. The inclusion of such radiation will intro- 
duce a pressure/density dependence via the ratio of the 
photoionization and recombination times. 

In a separate paper we will present computations of 
the integrated metal ion "cooling columns" produced in 
steady flows of cooling gas (e.g. Edgar & Chevalier 1986; 
Dopita & Sutherland 1996; Heckman et al. 2002) includ- 
ing the effects of "upstream" cooling radiation on the 
"downstream" ionization states. Such flows occur, for 
example, in the radiative (isobaric) cooling layers behind 
steady shock waves (Shull & McKee 1979; Draine & Mc- 
kee 1993). 

The outline of our paper is as follows. In §2 we write 
down the basic equations that we solve in our compu- 
tations, and we describe our numerical method. In §3 
we discuss isochoric versus isobaric cooling, and we de- 
rive critical cloud column densities and temperatures at 
which the dynamical transition from isobaric to isochoric 
cooling occurs. In §4 we present our computations of 
the ionization states as functions of gas temperature, 
for both CIE and non-equilibrium cooling. For time- 
dependent ionization we present results for both isochoric 
and isobaric evolutions. In §5 we present our computa- 
tions of the radiative cooling efficiencies. We compare 
the CIE and non-equilibrium cooling efficiencies, and de- 
termine the temperatures below which non-equilibrium 
effects become important for each assumed metallicity. 
In this paper we have generated a large number of fig- 
ures and tables containing results for metallicities rang- 
ing from 10~ 3 to 2 times solar. Some of the figures and 
tables appear in the main text, and some are available as 
online data files. Our data sets can be used to construct 
ion density-ratio diagnostics for radiatively cooling gas 
for a wide range of conditions. As an example, in §6 we 
discuss the evolution of the density ratios C IV/O VI ver- 
sus N V/O VI in radiatively cooling gas. We summarize 
in §7. 

2. BASIC EQUATIONS AND PROCESSES 

We are interested in studying the evolving ionization 
states of radiatively cooling gas, which is initially heated 
to some high temperature To > 5 x 10 6 K, and then 
cools from an initial (equilibrium) ionization state. In 
the absence of continued heating and ionization the ions 
recombine as the gas cools, and the overall ionization 
state decreases with time. If the gas cools faster than it 
recombines, non-equilibrium effects become significant, 
and the gas remains over-ionized throughout the cooling. 
Non-equilibrium ionization leads to a suppression of the 
cooling rates. We compute this coupled time-dependent 
evolution for clouds cooling at constant density or con- 
stant pressure, for a wide range of metallicities. Non- 
equilibrium effects are most important at high metallic- 
ity, where the cooling times are shortest. 

2.1. Ionization 
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In our numerical computations we follow the time- 
dependent ionization of an optically thin "parcel" of 
cooling gas, given an initial ionization state at an ini- 
tial temperature To. We consider all ionization stages 
of the elements H, He, C, N, O, Ne, Mg, Si, S, and Fe. 
The temperature-dependent ionization and recombina- 
tion processes that we include are collisional ionization 
by thermal electrons (Voronov 1997), radiative recombi- 
nation (Aldrovandi & Pequignot 1973; Shull & van Steen- 
berg 1982; Landini & Monsignori Fossi 1990; Landini & 
Fossi 1991; Pequignot, Petitjean, & Boisson 1991; Ar- 
naud & Raymond 1992; Verner et al. 1996), dielectronic 
recombination (Aldrovandi & Pequignot 1973; Arnaud & 
Raymond 1992; Badnell et al. 2003, Badnell 2006; Col- 
gan et al. 2003, Colgan, Pindzola, & Badnell 2004, 2005; 
Zatsarinny et al. 2003, 2004a, 2004b, 2005a, 2005b, 2006; 
Altun et al. 2004, 2005, 2006; Mitnik & Badnell 2004), 
and neutralization and ionization by charge transfer re- 
actions with hydrogen and helium atoms and ions (King- 
don & Ferland fits 1 , based on Kingdon & Ferland 1996, 
Ferland et al. 1997, Clarke et al. 1998, Stancil et al. 1998; 
Arnaud & Rothenflug 1985). Our atomic data set is 
the up-to-date compilation used in the steady-state pho- 
toionization codes ION (Netzer et al. 2005; and private 
communication) and Cloudy (Ferland et al. 1998). 

We assume that the gas is dust-free (see below) and ex- 
clude neutralization processes by dust grains. Our code 
can also account for photoionization by a steady or time- 
varying external radiation field. However, in the compu- 
tations presented here we exclude such radiation. 

The time-dependent equations for the ion abundance 
fractions, Xi, of element m in ionization stage i are, 

dxi/dt = Xi-i [qi-in c + IVi + fcj a i _ 1 n H + 

+ fc T*-l™He+] 

+x i+ i [a i+1 n c + kf i+1 n H a + fc^ 1 n He o] (1) 
-Xi [(qi + ai)n c + T, t + fc^nno + k^n^o 
+kf i n R + + kf?n He +] . 

In this expression, q-i and on are the rate coefficients for 
collisional ionization and recombination (radiative plus 
dielectronic), and k^, k^, k^, and k^ are the rate co- 
efficients for charge transfer reactions with hydrogen and 
helium that lead to ionization or neutralization. The 
quantities n B o , n H + , ri Hc o , n Bc + > an d n c are the particle 
densities (cm -3 ) for neutral hydrogen, ionized hydrogen, 
neutral helium, singly ionized helium, and electrons, re- 
spectively. Ti are the photoionization rates of ions i, due 
to externally incident radiation. Here, we set Ti = 0. 

For each element m, the ion fractions Xi = 
ni,m/(n-nA m ) must at all times satisfy 

where ni <m is the density (cm -3 ) of ions in ionization 
stage i of element m, tih is total density of hydrogen 
nuclei, and A m is the abundance of element m relative 
to hydrogen. The sum is over all ionization stages of the 
element. 

In our time-dependent code we can explicitly account 
for the effects of trapping of hydrogen and helium re- 
combination radiation by switching appropriately from 

1 See: http: / / www-cfadc.phy.ornl.gov/astro/ps / data/cx/hydro-gen 



"case A" to "case B" recombination, depending on the 
total cloud column density and temperature. As we dis- 
cuss in the Appendix, we find that such trapping has 
generally only a small effect on the ionization states of 
the metals in the radiatively cooling gas. We therefore 
present results assuming optically thin case A conditions 
throughout. We also neglect reabsorption of line and con- 
tinuum radiation emitted by metals in the cooling gas. 
In the Appendix we provide computational estimates of 
the maximal column densities for which reabsorption of 
the cooling radiation may be ignored. 

2.2. Cooling 

The ionization equations are coupled to an energy 
equation for the time-dependent heating and cooling, and 
resulting temperature variation. Here we are interested 
in clouds undergoing pure radiative cooling with no heat 
sources. 

For an ideal gas, the gas pressure P — Nk-gT/V ', and 
the internal thermal energy U — 3/2 Nk B T, where T 
is the gas temperature, V is the total volume, N is the 
total number of particles in the system, and k B is the 
Boltzmann constant. If dQ is an amount of heat lost (or 
gained) by the thermal electron gas, then 

dQ = dU + PdV = f 5 + S J {Nk B dT + k B TdN). (3) 
where 

JO, isochoric cooling, dV = , 

~ I 1 , isobaric cooling, dP = ^ ' 

For pure radiative cooling, with no heating, dQ = 
—n e nH^(T,Xi, Z)Vdt where A(T, Xi,Z) is the electron 
cooling efficiency (erg s -1 cm 3 ), which depends on the 
gas temperature, the ionization state, and the total abun- 
dances of the heavy elements specified by the metallicity 
Z. The electron cooling efficiency includes the removal 
of electron kinetic energy via recombinations with ions, 
collisional ionizations, collisional excitations followed by 
prompt line emissions, and thermal bremsstrahlung. For 
the low density "coronal limit" that we consider, A is 
independent of the gas density or pressure. 

In equation @ we do not include the ionization po- 
tential energies as part of the total internal energy 
(e.g. Schmutzler & Tscharnuter 1993) since dQ refers to 
the heat lost by the thermal electrons. In our defini- 
tion of A, therefore, the ionization potential energy that 
is released as recombination radiation does not appear. 
Only the kinetic energy of the recombining electrons con- 
tributes to the cooling efficiency 2 . On the other hand, ki- 
netic energy removed via collisional ionization is included 
in our A. If ionization potential energy is considered as 
part of the total internal energy, then collisional ioniza- 
tion does not lead to a net energy loss, since the kinetic 
energy removed is merely stored as potential energy. Ei- 
ther way of accounting for the energy losses leads to the 
same temperature versus time relation T(t). 

It follows from equation (J3J) that the gas temperature 

2 The recombination radiation energy e r = cke + ei, where eke 
is the thermal kinetic energy of the recombining electron, and ej is 
the ionization energy of the recombined ion. In our A we include 
ratnd^cfcjhjpil 
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declines at a rate given by (e.g., Kafatos 1973), 

dT n c nnA(T,Xi, Z) T dN 
~dt ~ (3/2 + s)nk B N~dt 

where n = N/V is the total particle density of the gas. 
Because of the PdV work, the temperature declines more 
slowly for isobaric (s = 1) cooling compared to isochoric 
(s = 0) cooling. 

The second term in equation j£| reflects the relative 
change in the total number of particles in the system. 
At temperatures > 10 4 K, where the hydrogen and he- 
lium remain ionized, the number of particles remains ap- 
proximately constant and this term may be neglected. It 
plays a role at lower temperatures as the hydrogen re- 
combines (Kafatos 1973). For a primordial helium abun- 
dance, yljje = TiHe/rva — 1/12 (Ballantyne, Ferland, & 
Martin 2000), the temperature declines at a rate 



TABLE 1 

Solar Elemental Abundances 



dT 
dt 



nA 



4.34 x (3/2 + s)/c E 



(6) 



for gas in which the hydrogen and helium are fully ion- 
ized. 

Equation JSJ can also be expressed as (e.g., Shapiro & 
Moore 1976), 

3dP hPdp 

2H-2-pA = - ncn * K{T ^ Z) (?) 

where p = fin is the mass density of the gas, and [i is 
the mean mass per particle. For hydrogen and helium 
fully ionized, and a primordial helium abundance, /i = 
16/27 toh where mn is the proton mass. 

Equations (JIJ and (|5jl imply that in the absence 
of external radiation (r^ = 0) the derivatives of the 
ion fractions with respect to temperature, dxi/dT = 
(dxi/dT)(dT/dt), are independent of the gas density or 
pressure. Because collisional ionization, electron-ion re- 
combination, hydrogen and helium charge-transfer reac- 
tions, and the gas cooling processes are all "two-body" 
interactions, with rates per unit volume proportional to 
n 2 , the density dependence divides out. The solutions 
for the ion fractions as functions of the instantaneous 
gas temperature, Xi(T), are therefore independent of the 
assumed gas density or pressure. 

Given a set of non-equilibrium ion abundances, Xi(T), 
and an assumed metallicity Z, we use the cooling func- 
tions included in Cloudy (version 06.02, Ferland ct 
al. 1998) to calculate A(T, Xi,Z), We present results 
for Z ranging from 10~ 3 to 2 times solar. For Z = 1 
we adopt the elemental abundances for C, N, O, Mg, 
Si, S, and Fe reported by Asplund et al. (2005a) for the 
photosphere of the Sun. An exception is Ne, for which 
we adopt the "enhanced" abundance recommended by 
Drake & Testa (2005) 3 . We list our assumed solar abun- 
dances in Tabled In all computations we assume a pri- 
mordial helium abundance A Bc = 1/12, independent of 
Z. The elements we include are the dominant coolants 
over the entire temperature range we study. Cooling by 
other elements is negligible. 

3 An enhanced Ne abundance may reconcile the "new" Asplund 
photospheric abundances with solar models and helioseismology 
(Bahcall, Basu, & Serenelli 2005; Anita & Basu 2005; but see As- 
plund et al. 2005b; Schmelz et al. 2005). See, however, Ayres et 
al. (2006). 



Element 


Abundance 




(X/H) 


Carbon 


2.45 x 10" 4 


Nitrogen 


6.03 x 10" 5 


Oxygen 


4.57 x 10" 4 


Neon 


1.95 x 10" 4 


Magnesium 


3.39 x 10" 5 


Silicon 


3.24 x 10" 5 


Sulfur 


1.38 x 10" 5 


Iron 


2.82 x 10" 5 



Non-equilibrium effects become important when the 
cooling time 

_ T 4.34 x (3/2 + s)k B T 
c ~\dT/dt\~ nA(T, Xi ,Z) [) 

becomes short compared to the recombination time for 
an ion i, 

— Z ' tot ' (9) 



where a* ot is the total recombination coefficient at 
the temperature T. The ionization state remains 
close to "collisional ionization equilibrium" (CIE) when 
t c /t r .i » 1 for the most abundant ions at the given 
temperature. Non-equilibrium effects become significant 
when t c /t It i < 1. 

Because different ions have different recombination 
times, the recombination-lag is not similar for all ions. 
Ions with longer recombination times, in particular the 
helium-like ions, may be expected to persist over a wider 
range of temperatures when out of equilibrium. 

Crucially, the ratios t c /t Tt i are independent of the gas 
density or pressure, so that departures from CIE are in- 
dependent of the assumed gas density or pressure, for 
both isochoric or isobaric cooling. However, because t c 
is shorter for isochoric cooling, non-equilibrium effects 
may be expected to be somewhat more pronounced for 
isochoric cooling compared to isobaric cooling. 

The primary, and essentially only, "free parameter" in 
the problem is the metallicity Z of the gas, through its 
control of A and the cooling time. A higher Z leads to 
enhanced metal line cooling, and larger departures from 
CIE, and vice versa. 

In our computations we assume that the gas is dust- 
free, and we do not consider cooling via gas-grain colli- 
sions, a process that can dominate the total cooling at 
high temperatures if a high dust mass can be maintained 
(Ostriker & Silk 1973, Draine 1981). The neglect of dust 
is appropriate when the dust sputtering destruction time 
scale, t s , is shorter than the initial cooling rate, t c , so that 
any initial dust is destroyed before appreciable gas cool- 
ing occurs. Here we rely on the computations of Smith 
et al. (1996) who found that t s /t c < 1 for T > 3 x 10 6 
K (see their Fig. 7). In our calculations we consider gas 
cooling from initial temperatures greater than this, and 
we assume that the dust is instantaneously destroyed. 
The gas then evolves with constant gas-phase elemental 
abundances specified by the assumed initial metallicity 
Z. 
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2.3. Numerical Procedure 

For the atomic elements that we include, equations |T|I 
and J5J are a set of 103 coupled ordinary differential 
equations (ODEs). In our numerical procedure we ad- 
vance the isochoric solutions in small pressure steps 
AP = e P where e < 0.05, and P is the current pressure, 
associated with a temperature T{P). For isobaric con- 
ditions, we advance the solutions in small density steps 
Ap = e p, where p is the current mass density, associated 
with a temperature T(p). 

For any T we compute the total cooling rate by passing 
the current (non-equilibrium) ionization fractions Xi(T) 
to the Cloudy cooling functions. We then use equa- 
tion (7} to compute the time interval At associated with 
pressure change AP (or density change Ap). We then 
integrate equations over the time interval At using a 
Livermore ODE solver (Hindmarsh 1983). When inte- 
grating the ionization equations , we assume that over 
the time step At, the pressure (or density) evolves lin- 
early with time. In the integrations, the estimated local 
errors on the fractional abundances were controlled so as 
to be smaller than 10~ 6 , 10~ 5 , and 10~ 4 for hydrogen, 
helium, and metal-ions, respectively. 

The above procedure is repeated with the new values of 
the ionization fractions, Xi(T+ AT), down to a minimum 
temperature Ti ow . Here we set 71 ow = 10 4 K. We ver- 
ify convergence by re-running the computation at higher 
resolution (smaller e), and confirming that the resulting 
ionization and cooling rates as functions of time remain 
unaltered. 

3. ISOCHORIC VERSUS ISOBARIC COOLING 

A gas cloud will cool isochorically (i.e., at constant 
mass-density and volume) when the cooling time as de- 
fined by equation JSJ is short compared to the dynamical 
time 

t d =- (10) 



where D is the cloud diameter, and c s — yJk^T / u is the 
sound speed. Cooling is isobaric (i.e., occurs at constant 
pressure and decreasing volume) when the cooling time 
is long compared to the dynamical time. 

The dynamical evolution is therefore determined by 
the ratio 



tc 

— ~ 4.34 x 

t d 



(3/2 + s)(fc B T) 3/2 



C 2 7iAD 



(11) 



Cooling is isochoric when t c /td << 1 and is isobaric when 
t c /td >> 1- Equation Ijlll) shows that for a given cloud 
size and density the cooling is isobaric if the tempera- 
ture is sufficiently high. This is because T 3 / 2 /A gener- 
ally increases rapidly with increasing temperature (see 
§5). Conversely, for a given temperature and density, 
the cooling is isochoric for sufficiently large D. 

For isochoric cooling, D and n remain constant with 
time. For these conditions t c /td oc T 3 / 2 /A, which gen- 
erally decreases as the gas cools. Hence, cooling that is 
initially isochoric remains isochoric. 

4 This package integrates initial value problems for stiff or non- 
stiff systems, and switches automatically between stiff and non-stiff 
methods as necessary. See: http://www.netlib.org/odepack/ 



For isobaric cooling, a spherical cloud with initial di- 
ameter Dq , density no , and temperature To , will contract 
to a size D = D (T/T ) 1/3 with density n = n T /T af- 
ter cooling to a temperature T. Thus, as the cloud cools 
and contracts t c /td oc T 13 / 6 /A. This ratio decreases 
rapidly as the gas cools, so that a transition to isochoric 
cooling must eventually occur. 

The transition from isobaric to isochoric cooling occurs 
at a critical cloud size, D tr , and temperature, Tt r , at 
which t c — td- It follows from equation lllj l that the 
transition occurs at a critical column density 



D tI n - 10 fc^V 172 ^ 

40 x T 6 2 ' t ° r 4 pc cm~ 3 



(12) 



where n is the gas density. Given that non-equilibrium 
effects reduce the cooling efficiencies by only factors 
a few (see §5) compared to CIE cooling, we use the 
CIE cooling efficiencies, A cq , in equation (|12|) to esti- 
mate Dtr and T tr . The resulting estimates for the crit- 
ical sizes and temperatures are then independent of ini- 
tial conditions. For the numerical evaluation in equa- 
tion (12) we use the power-law approximation A^ =1 = 
2.3 x 10 -19 7 1-0 - 54 erg cm 3 s" 1 for equilibrium cooling 
at solar metallicity (see §5). 

For a cloud that initially cools isobarically at pressure 
Po, from an initial diameter Do and temperature To, it 
then follows that 



nl/3 



io 4 /2 /i~ 1/2 T< 



A(ttr) 



= 4 x 10 7 Ti 



2.7 pc (cm 3 K) 



(13) 



6,tr (106 K )!/3 



which provides an implicit relation for the temperature, 
T t r j at which the transition to isochoric cooling occurs 
(cf. Edgar & Chevalier 1986). The numerical evaluation 
is for the power-law cooling efficiency at solar metallicity. 
In Figure H we plot Dtr" versus Tt r (left panel), and 

1 /3 

Ttr versus the parameter xo = DoPo/T (right panel) 
as given by equations ljT21) and ljT3|) . We present results 
for metallicities Z ranging from 10~ 3 to 2 times solar, 
assuming the CIE cooling efficiencies we compute in §5. 

In Figure HJi, the cooling for each Z is isochoric for 
temperatures and cloud column densities to the left and 
above of the curves. The cooling is isobaric to the right 
and below the curves. The dotted line in Figure^Ji is the 
power-law approximation given in equation (|12fl . The ar- 
rows in Figure indicate the direction in which the gas 
evolves as it cools. For isochoric cooling the column Dn 
remains constant, and the gas cools along horizontal tra- 
jectories. For isobaric cooling Dn oc T -2 / 3 , as indicated 
by the inclined arrows, showing that the gas must eventu- 
ally make the transition to isochoric cooling. For exam- 
ple, a condensation with diameter D = 10 kpc, density 
n = 10~ 5 cm 3 , and temperature T = 10 6 K (correspond- 
ing to a thermal pressure P = 10 cm -3 K) cools isobari- 
cally, since Dn = 0.1 pc cm -3 , whereas for this temper- 
ature D tr n = 40 pc cm~ 3 for Z = 1, or 10 3 pc cm~ 3 for 

z<io- 2 . 

In Figured, Ttr is the temperature at which the tran- 
sition from isobaric to isochoric cooling occurs, given an 
initial state defined by the parameter xq. In computing 
T tr we again use the equilibrium cooling efficiencies for 
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Fig. 1. — Critical column densities, Dun, and temperatures, Ttr, for the transition from isobaric to isochoric cooling, for clouds 
with met allic itics Z ranging from 10 — 3 to 2 times solar, (a) Transition column, Dun, versus transition temperature, Tt r , as given by 
equation 1121 . The dotted line is the power-law approximation for Z = 1. In the region above the curves the cooling is isochoric. Below 
the curves the cooling is isobaric. The arrows indicate the direction in which the gas evolves as it cools, (b) Transition temperature, Tt r , 

from isobaric to isochoric cooling, versus the initial state parameter xq = DqPq/Tq^ 3 , as given by equation 1131 . The dotted line is the 
power-law approximation for Z = 1. 



the different values of Z. The dotted line in Figure ^ 
shows the power-law approximation for Z = 1 as given in 
equation (T3J. For our above example (with D — 10 kpc, 
P = 10 cm 3 K, T = 10 6 K), the parameter x Q = 10 5 , so 
that the initial isobaric condensation will begin cooling 
isochorically at T tr = 10 5 K for Z = 1, or at 4 x 10 4 K 
for Z < 10~ 2 . 

As we show in the Appendix, a cooling gas cloud re- 
mains optically thin up to column densities that are much 
greater than the critical transition columns in Figure ^ 
Expressions (|12l) and i|13|) . and the curves in Figure Q are 
therefore of broad applicability for determining the iso- 
baric versus isochoric dynamical evolution for radiatively 
cooling clouds. 

4. IONIZATION FRACTIONS 

We have carried out computations of the collisionally 
controlled ionization states of the elements H, He, C, N, 
O, Ne, Mg, Si, S, and Fe, as functions of gas tempera- 
ture for three sets of assumptions. First we assume CIE 
imposed at all T. Then we consider the non-equilibrium 
ionization states as functions of the time-dependent tem- 
perature for radiatively cooling gas, for constant pressure 
and constant density evolutions. Our results are dis- 
played in Figure [21 and listed in tabular form in Tables 
2-4, and in additional tables and figures that we describe 
below. 

The left-hand panels of Figure [21 show our results for 
the CIE ion fractions x° q (T). Each panel displays the 
ionization stages of a particular element. In our CIE 
computations we have integrated the time-dependent 
equations to the equilibrium states reached at late 
times when the formation and destruction rates for each 
ion are equal. We performed these time integrations at 
fixed temperatures T ranging from 10 4 to 10 8 K. For 
purposes of comparison with our non-equilibrium com- 



TABLE 2 
CIE Ion Fractions 



Temperature 
(K) 


H°/H 




H+/H 




He /He 


1.00 x 10 4 


9.99 x 10~ 


l 


9.30 x 10- 


4 


1.00 x 10° 


1.05 x 10 4 


9.97 x 10- 


i 


2.88 x 10" 


3 


1.00 x 10° 


1.10 x 10 4 


9.93 x 10- 


i 


6.98 x 10" 


3 


1.00 x 10° 


1.15 x 10 4 


9.84 x 10" 


i 


1.51 x 10" 


2 


1.00 x 10° 


1.20 x 10 4 


9.69 x 10- 


l 


3.05 x 10- 


2 


1.00 x 10° 



Note. — (1) The complete version of this table is in the elec- 
tronic edition of the Journal. The printed edition contains only a 
sample. (2) CIE ion fractions are given for Z = 1. For T < 10 4 K 
the electron density, and associated ion fractions, depend on Z. 



putations, in Figure [21 we display results for tempera- 
tures between 10 4 and 5 x 10 6 K only. Results for higher 
temperatures, up to 10 s K where nearly all of the ions 
become fully stripped, are given in Table 2. We verified 
that the CIE solutions we obtained by time-integration 
are identical to those found by solving the set of alge- 
braic equations for xi obtained by setting dxi/dt — in 
equation as appropriate for steady state 5 . 

For an ionized hydrogen gas, the CIE ion fractions 
are independent of the metallicity Z ', and are "univer- 
sal" functions of the gas temperature T , as determined 
by the microscopic steady-state balance between electron 
impact collisional ionization, electron-ion recombination, 
and charge transfer 6 . Our CIE solutions, and in partic- 

5 We solved the algebraic equations using both Cloudy (Fcrland 
et al. 1998) and ION (Netzer et al. 2004) assuming identical in- 
put atomic data sets. These two "photoionization codes" assume 
steady-state conditions. 

6 For T < 10 4 K, hydrogen becomes neutral and the electron 
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ular the shapes and peak positions of the x^ q (T) curves, 
are in excellent agreement with previous computations 
available in the literature (see §1). Our CIE ion frac- 
tions differ only slightly from the widely quoted results 
of Sutherland & Dopita (1993). For example, as found 
by Sutherland & Dopita for carbon (see their Fig. 3), 
x(C 3+ ) peaks at a value of 0.3 around a narrow temper- 
ature range close to 10 5 K, whereas x(C 4+ ) is close to 
1 for a broad temperature range from ~ 1.5 x 10 5 K to 
- 7 x 10 5 K. The widths of the x eq {T) curves reflect the 
magnitudes of the ionization potentials of the ions. For 
example, C 4+ , N 5+ , and 6+ persist for a wide range of 
temperatures due to the increase in thermal energy scale 
that is required to remove the inner K-shell electrons. 

The right-hand panels of Figure [5] display our results 
for non-equilibrium ionization in radiatively cooling gas 
for solar (Z = 1) composition. For time-dependent cool- 
ing the behavior depends on the initial conditions. Here 
we assume that the gas is cooling from an initially hot, 
T > 5 x 10 6 K, equilibrium state. For such high tem- 
peratures the cooling times are long, and there is gener- 
ally sufficient time for the gas to reach CIE before ap- 
preciable cooling occurs. Our results are also listed in 
"on-line" Tables 3 and 4. Results for other metallicities 
are presented in additional electronic on-line figures and 
tables. In each panel in Figure the black curves are 
for isochoric cooling, and the light grey curves are for iso- 
baric cooling. Departures from CIE become apparent for 
T < 10 6 K, for which cooling becomes rapid compared 
to recombination. 

When departures from equilibrium occur, the gas at 
any temperature remains over-ionized compared to CIE, 
as recombination lags behind cooling. The recombina- 
tion lag affects all elements including hydrogen and he- 
lium. For example, in CIE the hydrogen is half-ionized 
at T = 1.6 x 10 4 K, and is fully neutral (99.9%) at 10 4 K. 
However, for non-equilibrium (isochoric) cooling for Z = 
1 the hydrogen is only 9% neutral at 10 K. The metals 
are affected in a similar manner, with Xi(T) remaining 
overionized compared to CIE at low temperatures. The 
recombination lags are largest for ions with small recom- 
bination efficiencies. Thus, He-like ions, such as C 4+ , 
N 5+ , and 6+ , which have long recombination times, 
persist to much lower temperatures in radiatively cool- 
ing gas compared to CIE. For example, for solar composi- 
tion z(C 4+ ) > 0.2 down to T ~ 1.4 x 10 4 K for isochoric 
cooling, and 2.3 x 10 5 K for isobaric cooling, whereas 
in CIE x eq (C 4+ ) becomes vanishingly small for T below 
~ 8 x 10 4 K. Similarly, the non-equilibrium abundance of 
the widely observed ion 5+ remains greater than 10% 
of its peak value of 0.11 occurring at 3.0 x 10 5 K, down 
to temperatures of 2.0 x 10 4 K for isochoric cooling, and 
1.0 x 10 5 K for isobaric cooling. In CIE, 5+ vanishes 
below - 2 x 10 5 K. 

Figure [3 shows that at high temperatures, where CIE 
is attained, the isobaric and isochoric Xi(T) curves con- 
verge as they must. However, for lower temperatures, 
the recombination lags are greater for isochoric cooling 
compared to isobaric cooling. This is due to the shorter 
isochoric cooling times (see equation [S]). The isochoric 

density depends on the assumed metallicity. At these low temper- 
atures the associated CIE metal ion fractions then do vary with 
Z. 



and isobaric ion fractions differ most substantially for 
ions with the largest recombination lags. The differences 
grow with increasing metallicity as the cooling times be- 
come shorter. 

Another important feature of the time-dependent ion 
distributions is the "double-peaked" behavior of the 
Xi(T) curves seen for many ions, e.g. C ++ , C 3+ , Si 3+ , 
S 4+ , and S 5+ . This behavior reflects the temperature 
dependence of dielectronic versus radiative recombina- 
tion (Kafatos 1973). For such ions, dielectronic recom- 
bination is the dominant removal mechanism from the 
high-temperature peaks to the central minima. At lower 
temperatures dielectronic recombination weakens, and 
the ion fractions rise again due to continuing recombina- 
tion from the persisting higher ionization stages. As the 
temperature continues to decrease, radiative recombina- 
tion becomes increasingly effective and the ion fractions 
drop again, leading to the low-temperature peaks. The 
double-peaked behavior is a time-dependent effect, be- 
cause the low-temperature peaks occur at temperatures 
for which the ions no longer exist in CIE. 

Our non-equilibrium Xi(T) curves are in qualitative 
agreement with previous computations (Kafatos 1973, 
Shapiro & Moore 1976, Sutherland and Dopita 1993, 
Schmutzler & Tscharnuter 1993). However, there are 
some significant differences in detail which are due, at 
least in part, to differences in the assumed atomic pro- 
cesses and rate coefficients. For example, the early stud- 
ies by Kafatos (1973) and Shapiro & Moore (1976) did 
not include charge transfer reactions with hydrogen or 
helium, which affect the ion fractions for T < 10 5 K. 
Differences between our results and the more recent com- 
putations of Sutherland and Dopita (1993) and Schmut- 
zler & Tscharnuter (1993) are likely due to differences in 
dielectronic rate coefficients (which control the "double- 
peak" behavior) and in the assumed solar abundances. 
For example, Sutherland & Dopita adopted the Anders 
& Grevesse (1989) abundances with 0/H Q = 8.5 x 10~ 4 , 
whereas we adopt the more recent Asplund et al. (2005a) 
value of 4.6 x 10 -4 (see Table^). Different assumed abun- 
dances lead to altered cooling times and departures from 
CIE. 

Our results for the non-equilibrium ion fractions, 
Xi(T), for metallicities Z equal to 10 -3 , 10~ 2 , 10 -1 , and 
2, are presented in electronic "on-line" Figures 3-6, and 
Tables 5-12. Here we consider two examples of how the 
assumed metallicity affects the ion fractions. In Figure 
0we display the 5+ and ++ distributions for the dif- 
ferent values of Z in isochorically cooling gas. At high 
temperatures (> 5 x 10 6 K) the cooling times are long 
for all Z, and the 5+ and ++ abundances converge 
to the CIE distributions (thick grey curves). Departures 
from CIE occur at lower temperatures. The departures 
are largest for large Z where the cooling times are short- 
est. For example, for Z = 2 the 5+ distribution is very 
broad, and this ion persists down to ~ 10 4 K, remain- 
ing within a factor 3 of its peak abundance of 0.1 at 
T = 1.7 x 10 4 K. For smaller Z the cooling times are 
longer, and the deviations from CIE become smaller. 

For very low metallicities (Z < 10~ 2 ), the metals con- 
tribute negligibly to the gas cooling (see §5), and the 
cooling times are independent of Z (e.g. Boehringer &; 
Hensler 1989). In the "low-Z" limit, the ion-fraction 
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Fig. 2. — Ion fractions Xi = rt; im /rtHA m , versus gas temperature. Each row displays results for a different element. Left hand panels 
are for collisional ionization equilibrium (CIE). Right hand panels are for non-equilibrium isochoric and isobaric cooling for Z = 1 times 
solar metallicity gas. Dark curves are for isochoric cooling. Light curves are for isobaric cooling. Panels show H, He, C, and N. 



Time-Dependent Radiative Cooling 



9 




10 



Gnat & Sternberg 




Time-Dependent Radiative Cooling 



11 




Fig. 2. — Continued. Panels show time-dependent ion fractions for Si, S, and Fc. 
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TABLE 13 
CIE Cooling Efficiencies 



Temperature 


A(Z = W- 3 ) 


A(Z = 10~ 2 ) 


(K) 


(erg cm 3 s — 1 ) 


(erg cm 3 s _1 ) 


1.00 x 10 4 


4.70 x 10" 24 


4.98 X 10" 24 


1.05 x 10 4 


7.62 x 10" 24 


7.79 x 10" 24 


1.10 x 10 4 


1.23 x 10 -23 


1.25 x 10 -23 


1.15 x 10 4 


1.97 x 10~ 23 


1.98 x 10" 23 


1.20 x 10 4 


3.08 x 10" 23 


3.09 x 10" 23 



Note. — The complete version of this table is in the electronic 
edition of the Journal. The printed edition contains only a sample. 



curves therefore converge to specific "primordial" 
forms, which may or may not approach the CIE distribu- 
tions. This depends on whether the \ow-Z cooling time 
is long or short compared to the ion recombination time. 
For example, Figure [7^, shows that at low Z, the 5+ 
fractions approach CIE. However, Figure 03 shows that 
++ remains out of equilibrium in the \ow-Z limit. This 
is an example of an ion whose abundance distribution 
can never reach equilibrium, for any Z, in radiatively 
cooling gas. Such ions are generally those that appear 
at T < 2 x 10 5 K, where the \ow-Z cooling times remain 
sufficiently short due to efficient cooling by hydrogen and 
helium Lya emissions (see §5). 

5. COOLING EFFICIENCIES 

We have carried out computations of the cooling effi- 
ciencies A(T,Xi,Z) (erg cm 3 s" 1 ) assuming CIE at all 
temperatures, and for time-dependent radiative cooling 
for isobaric and isochoric evolutions. Our results are dis- 
played in Figure|Sl and listed in Tables 13-15. The upper 
panel of Figure |H1 shows our CIE cooling functions, A oq , 
for T between 10 4 and 10 8 K, and Z from 1(T 3 to 2. The 
lower panel displays the non-equilibrium cooling curves. 

In computing the CIE cooling efficiencies we use the 
equilibrium ion fractions x° q (T) computed in §4 as in- 
puts to the cooling functions. Our CIE results reproduce 
the standard "cosmic cooling curves" presented in many 
previous papers in the literature (see §1). Small differ- 
ences compared to previous calculations are expected due 
to updates in the input atomic data and in the assumed 
gas phase abundances of the heavy elements. 

For Z > 0.1, the radiative energy losses between 
10 5 and ~ 10 7 K are dominated by electron impact 
excitations of resonance line transitions of metal ions. 
Above ~ 10 7 K, metal line cooling becomes less effec- 
tive, and bremsstrahlung radiation dominates. The low- 
temperature peak at 2 x 10 4 K is mainly due to hydro- 
gen Lya cooling, a process that becomes less effective at 
higher temperatures where the neutral hydrogen fraction 
becomes small. Recombinations and forbidden-line tran- 
sitions are minor contributors to the gas cooling. The 
familiar peaks and features in the cooling curves are due 
to specific metal line coolants that dominate at different 
temperatures. The Lya maximum is followed by peaks 
at 1.0 x 10 5 , 3.0 x 10 5 , 5.0 x 10 5 , and 1.5 x 10 6 , due 
respectively to resonance line transitions of carbon, oxy- 
gen, neon, and iron ions. For Z = 1 we find a maximum 
equilibrium cooling efficiency of 4.6 x 10~ 22 erg cm 3 s _1 



at T = 2.3 x 10 5 K. 

For Z < 0.01, metal cooling becomes negligible, and 
the energy losses are dominated by hydrogen and helium 
only. H, He and He + line emissions dominate at low T, 
and electron bremsstrahlung due to scattering with H + 
and He ++ at high T. The cooling peak at 10 5 K that 
appears when Z becomes small, is due to He + Lya. 

Between 10 5 and 10 8 K, our Z — 1 cooling function 
is well fit by the power-law expression (cf. Kahn et 
al. 1976), 

Af q =1 = 2.3 x 10~ 19 T-°- 5i erg cm 3 s" 1 . (14) 

This approximation is accurate to within a factor of 2.4 
over this temperature range (10 5 < T < 10 8 K), and 
to within a factor of 1.7 for 10 5 < T < 6.3 x 10 7 K. 
For Z > 1, and between 10 5 and ~ 10 7 K, A is linearly 
proportional to Z. For this range of parameters it follows 
from equations (|14fl and (JSJ that the cooling time t c ~ 
0.22 Tg 1 - 54 / (n Z) Myr for isochoric cooling, and a factor 
1.6 times longer for isobaric cooling. 

For Z < 0.01, the cooling efficiency is independent of 
Z. In this low-Z limit, the two-piece power law 

A z = io- 3 _ / 1-7 x 10~ 16 T- 1 - 29 , 8 x 10 4 < T < 8 x 10 
CC J " \ 3.8 x 10~ 26 T 34 , 8 x 10 5 < T < 10 8 K 

, , ( 15 ) 
(in units of erg cm 3 s ) reproduces the CIE cooling 

efficiency between 8 x 10 and 10 8 K to within a factor 

of 1.35. 

Figure displays our results for the non-equilibrium 
cooling efficiencies. In these computations we use the 
non-equilibrium ion fractions, Xi(T), obtained by solving 
the coupled equations (1) and (5), as presented in §4, to 
compute A(T,Xi, Z). For each Z we assume sufficiently 
large initial temperatures such that ionization equilib- 
rium is established before significant cooling begins. De- 
partures from equilibrium then occur as the temperature 
drops and the cooling times become short. For each Z 
the filled circles in Figure |SJd indicate the points where 
the cooling curves begin to deviate by more than 5% from 
CIE cooling. This is also where the isobaric and isochoric 
cooling curves (indicated for each Z by the labels "IB" 
and "IC") bifurcate. For Z = 1, departures from equi- 
librium set in for T < 10 6 K. The deviations from CIE 
cooling begin at higher temperatures for higher Z, and 
lower temperatures for lower Z. For temperatures above 
the "departure points" the curves in the upper and lower 
panels of Figure 00 are identical. 

There are several important differences between the 
non-equilibrium and CIE cooling curves. First, the dis- 
tinct peaks and features that appear for CIE are smeared 
out in the non-equilibrium cooling curves. This is due to 
the broader ion distributions, Xi(T), that occur for non- 
equilibrium cooling. Individual ions then contribute to 
the cooling over a larger temperature range. For exam- 
ple, for non-equilibrium cooling at high Z ions such as 
++ and Ne ++ persist below 10 5 K, where they con- 
tribute significantly to the cooling in addition to Lya. 
The Lya cooling peak is itself broadened because a high 
electron fraction is maintained down to low tempera- 
tures. 

Second, the non-equilibrium cooling efficiencies are 
suppressed, by factors of 2 to 4, compared to CIE cooling. 
For example, for Z = 1 the maximum non-equilibrium 




cooling efficiency, occurring at T = 2.1 x 10 5 K, is reduced 
to 2.6 x 10 -22 erg cm 3 s _1 , a factor of 1.8 less than the 
maximal CIE efficiency. The suppression occurs because 
the gas remains "over-ionized" as it cools, and the den- 
sities of the specific coolants that are most effective at 
each temperature are reduced. The more highly ionized 
species that remain present at each temperature gener- 
ally have more energetic resonance line transitions, and 
these are less efficiently excited by the thermal electrons 
(McCray 1987). For example, for non-equilibrium cool- 
ing in the \ow-Z limit the absolute and relative heights 
of the hydrogen and helium Lya peaks are reduced com- 
pared to CIE cooling. This is because the H° and He + 
fractions remain smaller at any temperature compared 
to CIE. 

Third, because isochoric cooling is faster (see equation 
[8]) with correspondingly greater recombination lags, the 
suppressions in the cooling efficiencies are larger, and the 
isochoric curves fall below the isobaric curves. 

The differences we find between CIE and non- 
equilibrium cooling, are in good agreement with the re- 
sults found by Sutherland & Dopita (1993; see their 
Fig. 13) and Schmutzler & Tscharnuter (1993; see their 
Fig. 2.2). For Z = 1 the power-law expression 

A non- cq = 5.6 x 1(T 20 T-° 46 erg cm 3 s" 1 (16) 
provides a good fit to the non-equilibrium cooling effi- 
ciency between 10 5 and 10 s K, and is accurate to within 
a factor of 2.2. For 10 5 < T < 5.5 x 10 7 K, it is accurate 
to within a factor of 1.5. For \ow-Z, departures from CIE 
cooling occur only below 2 x 10 5 K, so that the fit given 
by equation (|15l) remains valid for radiatively cooling gas 
at low metallicity, for T > 2 x 10 5 K. 

The suppression of the cooling efficiencies compared to 
CIE implies longer cooling times. The actual recombi- 
nation lag is therefore smaller than would occur if the 
gas were able to cool at the faster CIE rates. The depar- 
tures from ionization equilibrium are thus stabilized by 
the reduction in the cooling rates. 



6. DIAGNOSTICS 

For a uniformly cooling gas cloud, the line-of-sight col- 
umn densities of ions i of element m may be expressed 
simply as N™ — A m Xi(T)NH, where Nh is the column 
density of hydrogen nuclei, A m is the abundance of ele- 
ment m (relative to H), and Xj(T) are the temperature- 
dependent ion fractions, as computed in 21 Specific col- 
umn density ratios 

7V« A n Xj(T) [ ' 

may therefore be used as diagnostic probes of the ion- 
ization state, temperature, and also metallicity in radia- 
tively cooling gas. For CIE, the ion fractions are inde- 
pendent of the metallicity Z, and the column density 
ratios depend only on the gas temperature (and the rel- 
ative elemental abundances A m /A n ). However, for non- 
equilibrium cooling the ion fractions also depend on Z , 
and hence so do the column density ratios. 

Diagnostic diagrams for radiatively cooling gas for dif- 
ferent metallicities, may be constructed using the com- 
putational data we presented in §4 (Figs. 2-6 and Tables 
2-12). As an example 7 , in Figure El we display "cooling 
trajectories" for Nc iv/No vi versus y/No vi, for 
T ranging from 5 x 10 6 to 10 4 K. The gas temperature 
is represented by color along the curves, from hot (red) 
to cool (blue). Panel (a) shows the behavior for CIE, 
for which the results are independent of Z. Panels (b)- 
(f) are for non-equilibrium isochoric cooling, for Z from 
10~ 3 to 2. The CIE trajectory is reproduced as the grey 
curves in panels (b)-(f). 

For the Z=l model in Figure ED we also display, as ob- 
servational examples, the data points and upper limits 
presented by Fox et al. (2005) for the multiphased high- 
velocity cloud absorbers towards the background sources 

7 Diagnostic diagrams for any combination of ion-ratios 
may be constructed automatically using our web tool at 
http: / /wise-obs. tau.ac.il/~orlyg/cooling/ 
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Fig. 8. — Cooling efficiency (erg cm s 1 ) versus temperature 
for Z from 10 — 3 to 2 times solar, (a) Equilibrium cooling. The 
dominant cooling elements at various temperatures are indicated 
near the curves. The dotted lines are our pow er-l aw f its for Z = 
1 and Z = 10 — 3 , as given by equations 1141 and 1151 (b) Non- 
equilibrium cooling. For each Z, two curves are displayed: the 
lower curve is for isochoric (IC) cooling, and the upper curve is 
isobaric (IB) coolin g. T he dotted line is our power-law fit for Z = 1 
as given by equation lltil The filled circles indicate the temperatures 
where time-dependent cooling first deviate from the equilibrium 
cooling by 5%. 



HE 0226-4110 and PG 0953+414 (black crosses). The 
plot also includes the Galactic halo average (blue cross) 
A c iv /N vi = 0.6±0.47 and A N Y /N V i = 0.12±0.07 
(Zsargo et al. 2003, as quoted by Fox et al.). We also 
show the regions demarcated by Fox et al. for the pre- 
dictions for a variety of ionization mechanisms, includ- 
ing turbulent mixing layers (TMLs; Slavin, Shull, & 
Begelman 1993), conductive interfaces (CIs; Borkowski, 
Balbus, & Fristrom 1990), cooling flows (CFs; Edgar 
& Chevalier 1986) and shock ionization (SI; Dopita & 
Sutherland 1996). 

To set the absolute scale of our predicted ion columns, 
in Figure ITU1 we plot Ac iv/Nq vi versus Ao vi/Ah,20 
for each Z , where Ah,2u is the hydrogen column density 
in units of 10 20 cm~ 2 . For the CIE calculation we assume 



Z = 1. At any temperature, the O VI column scales lin- 
early with the assumed Ah, whereas the column density 
ratios are independent of Ah- 

Of the three ions we are considering, 5+ is the most 
highly ionized, so that C IV/ O VI and N V/O VI are 
small at high temperature, and become large at low tem- 
perature. Figures ED and ^| show again that for time- 
dependent radiative cooling the gas remains more highly 
ionized compared to CIE. As Z becomes large and the re- 
combination lags grow, the ion ratios remain small even 
at low temperatures. The cooling trajectories are there- 
fore confined to a narrower range of ion ratios for non- 
equilibrium cooling. Furthermore, in our example, it 
is evident that for a given N V/O VI ratio, the associ- 
ated C IV/O VI ratio is larger for non-equilibrium cool- 
ing compared to CIE. For example, a cloud for which 
An v/Nq vi — 3 and Ac iv/No vi — 30 would be in- 
consistent with gas at CIE at any temperature, yet con- 
sistent with isochoric radiatively cooling gas at a tem- 
perature near ~ 1.3 x 10 5 K. It is clear that assuming 
CIE when interpreting such data may lead to significant 
inaccuracies. For example, for an isochorically cooling 
cloud with Z ~ 1, absorption line measurements yield- 
ing Ac iv /No vi ~ 1 and An v /No vi ~ 1 would be 
interpreted as indicating T ~ 2 x 10 5 K assuming CIE, 
whereas in fact these ratios occur at a much lower tem- 
perature ~ 10 4 K. 

Interestingly, many (but not all) of the Fox et al. (2005) 
data points lie close to our Z = 1 cooling trajectory (see 
Figure Efc) . At lower metallicities the column density 
ratios approach the CIE trajectory which lies well below 
the data points. However, as indicated by Figure|5|3, even 
for Z = 1 further constraints may be required to distin- 
guish radiatively cooling gas from ionization occurring in 
conductive interfaces or turbulent mixing layers. 

7. SUMMARY 

In this paper we present new computations of the equi- 
librium and non-equilibrium cooling efficiencies and ion- 
ization states for low-density radiatively cooling gas, con- 
taining cosmic abundances of the elements H, He, C, N, 
O, Ne, Mg, Si, S, and Fe. In these calculations we assume 
pure radiative cooling, with no gas heating or photoion- 
ization by external sources. We present results for gas 
temperatures, T, between 10 4 and 10 8 K. We assume that 
the gas is dust free, and we consider metallicities Z rang- 
ing from 10~ 3 to 2 times the elemental abundances in the 
Sun. We carry out our computations using up-to-date 
rate coefficients for all of the atomic recombination and 
ionization processes, and the energy loss mechanisms. 

For temperatures below ~ 5 x 10 6 K, where ion-electron 
recombination lags significantly behind the cooling, we 
explicitly solve the coupled time-dependent ionization 
and energy loss equations for the cooling gas. For such 
gas we assume that the cooling is from an initially hot 
equilibrium state. The basic equations and our numer- 
ical method are presented in §2. We calculate the non- 
equilibrium cooling efficiencies for constant pressure (iso- 
baric) and constant density (isochoric) evolutions. De- 
partures from collisional ionization equilibrium (CIE) are 
slightly smaller for isobaric cooling. 

In §3 we consider the conditions for isochoric versus 
isobaric cooling. We compute the critical column densi- 
ties and temperatures at which the cooling time becomes 
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Fig. 9. — Column density ratios Nq iv/No vi versus JVn v/^O VI for radiatively cooling gas. Gas temperature is indicated by color 
along the trajectories, from hot (red) to cool (blue). Temperature versus color legend is on the right, (a) displays the CIE trajectory, shown 
again as grey curves in the other panels, (b)-(f) are the trajectories non-equilibrium isochoric cooling for different values of the metallicity 
Z. The data points (Fox et al. 2005) show the ionic ratios observed in metal absorbers toward HE 0226-4110 and PG 0953+414 (black), and 
in the Galactic halo (blue). The boxes show model prediction, as demarcated by Fox et al., for turbulent mixing layers (TMLs), conductive 
interfaces (CIs), cooling flows (CFs) and shock ionization (SI). 



short compared to the dynamical time, and the transi- 
tion from isobaric to isochoric evolution occurs. These 
results are displayed in Figure ^ and are based on the 
cooling efficiencies we present in §5. 

Because we exclude photoionization by external radia- 
tion, both the equilibrium and time-dependent behavior 
is independent of the assumed density or pressure. The 
primary parameter for radiative cooling is the metallicity. 
Departures from equilibrium are largest for high Z where 
the ratios of the cooling and recombinations times are 
smallest. At very low metallicity (Z < 10 -2 ) the cooling 
rates, and hence also the recombination lags, are inde- 
pendent of the metallicity. In this limit any departures 
from CIE approach a universal "primordial" form. The 
results of our equilibrium ionization calculations are dis- 
played in Figure 2 in §4. This figure also shows the ion- 
ization states for non-equilibrium isobaric and isochoric 
cooling, for Z = 1. Results for other metallicities are pre- 
sented in on-line tables and figures, as summarized in Ta- 
ble 16. Our computational data set can also be found at 



our website http://wise-obs.tau.ac.il/~orlyg/cooling/ 

Figure [8] in §5 displays our results for the equilibrium 
and non-equilibrium cooling efficiencies for Z between 
10~ 3 and 2. We provide simple, or two-piece, power-law 
fits (equations [14], [15], and [16]) for the equilibrium 
and non-equilibrium cooling efficiencies at high and low 
Z . Non-equilibrium cooling is suppressed, by factors of 
2 to 4, relative to CIE, because of the generally higher 
ionization state of the rapidly cooling gas. The famil- 
iar peaks in the CIE cooling curve are smeared out for 
non-equilibrium cooling, because individual ions persist 
for a broader range of temperatures. Overall, our results 
are in good agreement with previous such calculations. 
Detailed difference are mainly due to differences in the 
input atomic data, and assumed abundances of the heavy 
elements. We make some explicit comparisons with pre- 
vious computations in §4 and §5. 

Ion ratios are useful as diagnostic probes. In §6 we 
discuss one example, N V/O VI versus C IV/OVI, and 
show how this ratio evolves in radiatively cooling gas, and 
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TABLE 16 
Ionization and Cooling Data 
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is optically thin. In the Appendix, we provide numeri- 
cal estimates for the maximal cloud column densities for 
which this assumption remains valid. We also investigate 
how reabsorption of hydrogen and helium recombination 
radiation in optically thick clouds alters the cooling rates 
and associated ionization states in non-equilibrium cool- 
ing gas. For high metallicity the effects of trapping are 
small, but become more significant (~ factors of 2) for 
low Z gas. 

In a companion paper we will present computations 
of the metal ion "cooling columns" produced in steady 
flows of radiative cooling gas, such as occur in post shock 
cooling layers, including the effects of "upstream" cool- 
ing radiation on the "downstream" ionization states. We 
will also consider how photoionization by external back- 
ground radiation fields alters the ionization and thermal 
evolution of radiatively cooling gas such as we have con- 
sidered here. 



how it can be used as a probe of metallicity for realistic 
non-equilibrium conditions. 
In our computations we assume that the cooling gas 
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APPENDIX 

REABSORPTION OF DIFFUSE RADIATION FIELDS 

In the computations presented in §4 and §5 we have assumed that the cooling gas is optically thin, and that 
reabsorption of line and continuum radiation emitted by the cooling gas may be ignored. In this Appendix we provide 
numerical estimates for the maximal cloud column densities for which the optically thin assumption is justified. We 
then examine how trapping of hydrogen and helium recombination radiation in optically thick clouds affects the 
time-dependent cooling efficiencies and ionization states. 

For sufficiently large cloud column densities, absorption of the internally generated "diffuse radiation" can alter the 
ionization state in such a way as to enhance or reduce the net cooling rate. To estimate the critical column densities 
at which the cooling rates are altered, we use Cloudy to compute the local CIE cooling rates as functions of the total 
hydrogen column density, JVh, for one-dimensional constant temperature slabs. We set the external radiation field 
to zero, so that the radiative transfer handled by Cloudy is for the internally generated diffuse fields only. For each 
temperature, T, and assumed metallicity Z, we identify the critical column density, iVjj C rit (T, Z) , at which the local 
cooling efficiency at the cloud center first deviates by 50% from the optically thin cooling rate at the cloud edge. 

The critical column densities defined in this way are sensitive to the gas temperature and associated ionization 
states. Here we are no t inte rested in the detailed fluctuations of -ZVncrit with gas temperature, but rather with the 
broad trends. In Figure lAl II we display a smoothed representation of iVjjcrit versus T. The solid curve, appearing for 
T < 6 x 10 4 K is for the entire metallicity range 10~ 3 < Z < 2. For higher temperatures, the behavior is sensitive 
to Z, and we display th e crit ical columns by dashed lines for selected metallicity ranges. For temperatures where no 
curves appear in Figure lAl II the critical columns exceed 10 24 cm -2 . 

Three distinct regimes appear in Figure lAllI For 10 4 < T < 6 x 10 4 K, a transition from "case A" to "case B" 
hydrogen recombination occurs as the cloud becomes optically thick. Photoionization of hydrogen by H + and He + 
recombination radiation then alters the dominating Lya cooling rate. The Lya cooling rate is proportional to the 
product x c (l — x c ), where x D is the electron fraction. For T < 1.5 X 10 4 K, the neutral fraction (1 — x e ) « 1, whereas 
x e increases in the transition to case B. This leads to increased Lya cooling. At these low temperatures we find 
that ATjjcrit ~ 10 17 cm" 2 . At higher temperatures, T > 2.5 x 10 4 K, the hydrogen is largely ionized and x e » 1 
so the transition to case B mainly reduces the neutral fraction (1 — x e ). This leads to decreased Lya cooling. For 
T = 2.5 x 10 4 K, we find iVncrit ~ 5 x 10 19 cm~ 2 . Between these two temperatures, x e (l — x e ) remains approximately 
constant in the transition to case B, and a local maximum of iVncrit = 4 x 10 20 cm -2 appears at T w 2 x 10 4 K. Above 
T~6x 10 4 K A^Hcrit rises sharply because the hydrogen becomes so highly ionized that Lya cooling is no longer the 
dominant coolant. 

Similarly, absorption of Hc ++ recombination radiation by He + ions alters the He + Lya cooling rate at temperatures 
near 10 5 K where helium makes the transition from singly to doubly ionized form. However, this effect is only important 
for low metallicity (Z < 10~ 2 ) clouds, where He + Lya is a significant coolant. For such clouds, iVncrit = 1-5 x 10 22 cm~ 2 
at T = 1.6 x 10 5 K. 

Finally, for 10 6 < T < 5 x 10 6 K, we find that 7V Hc rit ~ 10 21 cm -2 in high metallicity (Z > 1) clouds. At these 
temperatures, photoionization by numerous energetic metal emission lines ionize 7+ , and shift the peak of the iron 
ion distribution from Fe 11+ -Fe 14+ to Fe 13+ -Fe 16+ . This reduces the cooling rate, because 7+ , Fe 11+ , and Fe 12+ 
are dominant coolants at these temperatures. For lower bremsstrahlung cooling plays a more dominant role, and 
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photoionization of the oxygen and iron i ons is less important, so Njicnt is much larger. 

Our computations for -/Vkcrit in Figure lAl II are for CIE conditions. We now examine how trapping of hydrogen and 
helium recombination radiation alters the cooling efficiencies and non-equilibrium ionization states for time-dependent 
cooling. We consider a series of isochoric model clouds with total column densities ranging up to iVn = 10 24 cm~ 2 . 
We assume that the H + , He + , and Hc ++ recombination photons are reabsorbed on-the-spot when the photoionization 
optical depths at the H° , He and He + ionization thresholds exceed unity respectively. The critical temperatures below 
which the clouds become optically thick are higher for clouds with larger total column densities, since lower neutral 
hydrogen and helium fractions are required. When the clouds become optically thick we switch from case A to case B 
recombination in computing the evolution of the hydrogen and helium ionization states 8 . 

In Figure IA12I we plot the neutral hydrogen fractions, xn(T), and cooling efficiencies, A(T), for case A and case B 
recombination. We display results for Z = 1 (left-hand panels), and for Z = 10~ 3 (right-hand panels). To illustrate 
the effects on the metal ion distributions we also plot the C + , C ++ , and C 3+ fractions. 

The neutral hydrogen fractions are, of course, reduced in optically thick clouds. The vertical lines connecting the 
case A and case B curves indicate the temperatures at which, for a given total cloud column TVh, the cloud becomes 
optically thick to hydrogen recombination radiation. For example, for iVn = 10 cm , this occurs at T = 5 x 10 4 K. 

For Z = 1 the differences between case A and case B cooling are small, and are only apparent at low temperatures 
where Lya cooling dominates. For non-equilibrium cooling the electron fraction remains large down to 10 4 K. Therefore, 
down to this temperature Lya cooling is reduced by the reduction in the neutral H fraction. (This is opposed to the 
situation in CIE where, as discussed above, for 10 4 K gas the Lya cooling is enhanced by an increased electron 
density in optically thick clouds.) For Z = 1, the slightly altered cooling efficiencies in optically thick clouds lead to 
only very small changes in the non-equilibrium metal ion abundances, as is illustrated for the carbon ions in Figure 
IA12I Our conclusions for Z = 1 are consistent with the results of Kafatos (1973) and Shapiro & Moore (1976), 
who considered the transition from optically thin to thick (for hydrogen recombination radiation only) at a single 
temperature T = 3.5 x 10 4 K (corresponding to N-r ~ 1.3 x 10 20 cm -2 ). 

For Z = 10~ 3 , hydrogen and helium are already the dominant emission line coolants below ~ 10 6 K, and the 
transition to o ptically thick conditions has a larger effect. The cooling efficiencies are reduced by up to a factor ~ 2. 
In Figure EH we draw three cooling curves for Z = 10 3 . The upper curve is for optically thin case A cooling. The 
middle curve (labeled "He + case B" ) is for clouds that are still optically thin to H + and He + recombination radiation, 
but are thick to He ++ recombination photons. The lower curve is for case B clouds that are thick to the H + and 
He + recombination photons as well. The vertical lines connecting the curves indicate the transition temperatures for 
different cloud columns N-r. For example, a 10 21 cm' 2 cloud shifts from case A to "He + case B" at ~ 2 x 10 5 K, and 
then becomes thick to H + and He + recombination radiation ("H/He case B") at ~ 5 x 10 4 K. 

For \ow-Z gas the overall cooling rates are small, and hence so are the departures from CIE. Therefore, even though 
the cooling efficiencies are significantly altered by trapping of the recombination radiation, the resulting effect on the 
ion fractions is not very large. As shown for C + , C ++ , and C 3+ the effects are most significant at low temperatures 
where the departures from CIE are largest. 
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